Mixed state of a dirty two-band superconductor: application to MgB 2 
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We investigate the vortex state in a two-band superconductor with strong intraband and weak 
interband electronic scattering rates. Coupled Usadel equations are solved numerically and the dis- 
tributions of the pair potentials and local densities of states are calculated for two bands at different 
values of magnetic fields. The existence of two distinct length scales corresponding to different 
bands is demonstrated. The results provide qualitative interpretation of recent STM experiments 
on vortex structure imaging in MgB2. 



A very peculiar feature of the recently discovered su- 
perconductor MgB2 0| is the multigap nature of the 
superconducting state. The possibility of such a state 
was first predicted in ||, || for a multiband superconduc- 
tor with large disparity of the electron-phonon interac- 
tion for the different Fermi-surface sheets. Various as- 
pects of multiband superconductivity, in particular the 
role of impurity scattering, were discussed theoretically 
in j| |, |, 0. For MgB 2 , the two-band model was first 
suggested in ||, ||] . On the basis of first-principles calcula- 
tions of the electronic structure and the electron-phonon 
interaction, it was argued that superconductivity in this 
compound resides in two groups of bands: the group of 
two strongly superconducting er-bands and the group of 
two weakly superconducting 7r-bands. Quantitative pre- 
dictions for T c , energy gaps, specific heat (h], |ll| and 
tunneling [Q were made recently for MgB2. 

Signature of two energy gaps was observed in Nb 
doped SrTi0 3 [[b|. But to date, only in MgB 2 exis- 
tence of two distinct gaps has been most clearly demon- 
strated. A large number of experimental data, in partic- 
ular tunneling jL4| [l5| and point contact measurements 
|l6| , |l7| , |Ls| , [lfjf and heat capacity measurements [^0| , di- 
rectly support the concept of a double gap MgB2. It was 
argued in Ref. |2l) that the unexpectedly weak correla- 
tion between T c and the resistivity can be reconciled with 
the two-band model, if one assumes that the interband 
impurity scattering remains weak even in samples with 
strong intraband impurity scattering in the 7r-band. 

Two-band superconductivity in MgB2 offers new in- 
teresting physics. For example, it was demonstrated 
that the anisotropics of the upper critical fields and the 
London penetration depths are different and have oppo- 
site temperature dependencies ^2|. Recently, the c-axis 
Abrikosov vortex structure in MgB2 was studied by STM 
p3| . Several important observations have been made: 
large vortex core size compared to estimates based on 
H C 2, the absence of zero-bias singularity in the core and 
the rapid suppression of the apparent tunneling gap by 
magnetic fields much smaller than H C 2 . Important prop- 
erty, that is essential for understanding these findings, is 



that c-axis tunneling in MgB2 probes mainly the weakly 
superconducting 7r-band [ jl5| . 

In this paper, we provide a quantitative model for 
the vortex structure in a two-band superconductor. We 
demonstrate the existence of two different spatial and 
magnetic field scales, consistent with the data in f23fl . 

We consider a two-band superconductor with weak in- 
terband impurity scattering and rather strong intraband 
scattering rates exceeding the corresponding energy gaps 
(dirty limit). In this case the quasiclassical Usadel equa- 
tions are applicable within each band. The vortex 
structure in single-band dirty superconductors was stud- 
ied extensively in the framework of the Usadel equations 
p5[ pfj[ . To describe the mixed state in the considered 
case, one can generalize the approach [^5[ |26| and write 
down the system of coupled Usadel equations 
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where a = 1,2 is the band index, A is the matrix of 
effective coupling constants (to be defined below), T> a 
arc diffusion constants, which determine the coherence 
lengths £ Q = ^/V\J2ttTc, G a and F a are normal and 
anomalous Green's functions connected by normalization 
condition G 2 a + F*F a = 1, A a is the pair potential and 
u) = (2n + 1)ttT are Matsubara frequencies. Bearing in 
mind the application to MgB2, in our notations index 1 
corresponds to er-bands and index 2 to 7r-bands. 

Note that in the considered case of weak interband 
scattering the Green's functions in different bands are 
coupled only indirectly, via the self-consistency equation 
(||). As will be shown below, this fact leads to the ex- 
istence of two different length scales in different bands 
and, as a consequence, two magnetic field scales appear 
which are directly accessible experimentally. Such a sit- 
uation has never existed in the field of vortex physics. 
This is in contrast to the usual proximity effect in real 
space (e.g. N/S multilayers), where different energy and 
length scales exist in spatially separated N, S layers. 



Let us study the case when magnetic field is oriented 
along c-axis Further, we neglect in-plane anisotropy and 
adopt a circular cell approximation for the vortex unit 
cell J2J| (see inset in Fig. [I] ). We also assume a large 
Ginzburg-Landau parameter k, k ^> 1 (this assumption 
is fulfilled in MgB 2 ) and consider magnetic fields signifi- 
cantly larger than the lower critical field, so that we can 
neglect variations of the magnetic field. To facilitate the 
analysis, we introduce reduced variables: we will use irT c 
as a unit of energy, and £1 = ^T>i/2ttT c as a unit of 
length. The distribution of superfluid momentum within 
the unit cell of the lattice is then given by 

p = 1/r - r/r 2 c , r\ = H x /H, H a = T c $ /V a (3) 

where r is the distance from the center of a vortex core. 

Using 9— parametrization (F a — smO, G a — cos9 a ) 
the Usadel equations and the self-consistency conditions 
can be rewritten in the form 
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with k\ = 1, fcf = Vx/V 2 , ui = t(2n + l),t = T/T c . The 
matrix W a p is related to the coupling constants A Qj a as 
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Wi = — , VF 2 = — 
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W12 = Aia/Det, TF 2 i = A 2 i/Det 
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where A = (An - A 22 )/2, Det = A U A 22 - Ai 2 A 2 i @. 
Note that only three constants are independent since 
W1W2 = W12W21. 

Partial local densities of states (DoS) N a (e,r), which 
are accessible in tunneling experiments, can be obtained 
from Q a (u>, r) using analytic continuation 

N a (s, r) = Re [cos Q a (iu) — ► e + iS, r)] (7) 

The above set of equations fully defines vortex 

core structure in a diffusive two-band superconductor. In 
general, numerical solution is required to determine the 
behavior of the pair potentials and DoS in both bands. 
The problem simplifies near the upper critical field when 
Eqs. (^) can be linearized 

d 2 9 a + ±-d r 9 a -(l/r - r/r 2 f 9 a -k 2 a Lu8 a = ~k 2 a A a . (8) 

These equations have exact solution |25| 
A a = Ao Q rexp (-r 2 /2r 2 ) , 9 a = 9 0a rexp (-r 2 /2r 2 c ) , 
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FIG. 1: Local density of state for 7r-band at different points 
of vortex lattice unit cell at h = 0.1 for T>\ = O2 = L Inset 
illustrates the circular cell approximation and shows points at 
which the spectra are taken. 



giving the relation 
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Substituting this result into the self-consistency equa- 
tions we obtain 
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-W 2 i A01 + W 2 A 02 = A 02 fin - t -g (J^j \ (10b) 

where g(x) = ^(1/2 + a;) — ^(1/2) and ip( x ) 1S a digamma 
function. This gives the equation for H c2 
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and relation between A01 and Ao 2 near H C 2 
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In the single band case, it follows directly from Eq.(lOa) 
that the upper critical field H s c2 is given by the standard 
Maki - de Gennes equation 



Hl/t) = g[H s c2 /{tH 1 ) 



(13) 



Application to MgB 2 . The electron-phonon inter- 
action in MgB 2 was calculated from first principles in 
In Ref. (H the matrices of the electron- 
phonon coupling constants A<y and the renormalized 
Coulomb pseudopotentials fi*^ were derived for the ef- 
fective two-band model. In this paper we will use these 
results and define the effective constants Ay = A y — fj,*j in 
the weak coupling model, neglecting the strong-coupling 
corrections, which is a reasonable approximation for our 
purpose. The corresponding numerical values are [ fLl| 

A u w0.81, A 22 « 0.278, A 12 « 0.115, A 21 « 0.091, (14) 



FIG. 2: Spatial dependencies of pair potentials ((a) and (b)) 
and partial DoS at E = ((c) and (d)) for isolated vortex for 
two ratios D x jD 2 : £>i = 0.2D 2 and D 1 = D 2 



FIG. 3: Field dependencies of maximum pair potentials ((a) 
and (b)) and averaged DoS at e = ((c) and (d)) for two 
ratios D 1 /D 2 : D x = 0.2D 2 and D y = D 2 



from which we obtain values of W a p used in numerical 
calculations, 

Wi ~ 0.088, W 2 ~ 2.56, W 12 ~ 0.535, W 2 i « 0.424. (15) 

With fixed coupling constants the overall behavior is de- 
termined by the ratio of the diffusion constants T>i/T> 2 , 
which is not known at present and may depend on the 
type of scatterers. Available estimates of scattering rates 
p9[ suggest that typically T>i > Z? 2 - However we found 
that to describe the experimental vortex core structure 
we have to take T>i < T> 2 (as will be discussed below, this 
apparent contradiction can be resolved by going to £i/£ 2 
ratio). Therefore, we present calculations for two values 
of the ratio V x jV 2 : V x jV 2 = 0.2 and V x jV 2 = 1. 

The magnetic field is measured with respect to the 
single-band upper critical field of the cr-band, h = 
H/H^ 2 (t), where H s c2 is given by equation (|l3|). For 
the case W\ <C W 2 realized in MgB 2 , the upper criti- 
cal field is mainly determined by the strong band. Small 
correction due to the weak band can be found from Eq. 
|lT| using expansion with respect to the small parameter 
5i2 = W\/W 2 . In particular we found very simple ex- 
pressions for the slope of H c2 at T c and H c2 (0): 
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H c2 (0) « ff c s 2 (0) (1 + S 12 In {H 2 /Hx)) . 

With the above parameters, we numerically solved 
equations (ji|) and ( |5a| ) for different magnetic fields. Fig. 
|l| shows an example of local DoS for the 7r-band at dif- 
ferent points of the vortex unit cell. One can see that in 
the center of the core there is no zero-energy peak in the 
DoS in the core usually observed in clean superconduc- 
tors p7| . This property is a consequence of the dirty limit 
in the tt— band. As one can expect, the most pronounced 
dependence on energy is observed at the boundary of the 
vortex unit cell (curve 3 in Fig. [J) . One can see that the 
DoS is peaked at an energy about 3 times smaller than 
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FIG. 4: Ratio of the apparent vortex size as defined in 
Fig. ^c to the coherence length £ C 2 extracted from the upper 
critical field, £ c2 = \J <E>o / (2ttH c2 ) plotted vs ratio of the typ- 
ical len gth scales in two bands, £2/^1 = \J"Di /Th. According 
to Ref. [E3[ for MgE>2 single crystal the ratio £ v /£c 2 is around 
3. 

^maxi- This peak corresponds to the small energy gap 
in the second band. 

We study structure of an isolated vortex by solving 
the Usadel equations at very small field (h = 0.002). Fig. 
||a,b shows the spatial dependence of the pair potentials, 
Ai(r) and A 2 (r), and their ratio for i = 0.1 for two cases: 
£>i = 0.2D 2 and D x = D 2 . Fig. §c,d shows the DoS at 
zero energy, Ni(Q,r) and N 2 (0,r). One can see that in 
the case of D\ = 0.2D 2 the pair potential and the DoS in 
the 7r-band demonstrate qualitatively different behavior. 
The pair potentials approach their bulk values A Qj o a t 
the length scale set by the strong band. Ai reaches half 
of A^o at r = 2.15£i and A 2 reaches half of A 2i o at a 
somewhat larger length scale, r = 3.44£i. 7r-band DoS, 
N 2 (0, r), has significantly longer range: it drops to 0.5 at 
r = 6.35£i. Though the above numbers correspond to the 
specific choice of parameters for the coupling matrix Ay , 
the large core size in the weakly superconducting band is 
the general property of a two-band superconductor. 

The two typical sizes of the isolated vortex determine 
the two typical field scales. Fig. |^ shows the field depen- 
dence of the maximum values of the pair potentials at 
the boundary of vortex unit cell ((a) and (b)) and DoS 
at e = averaged over the unit cell ((c) and (d)). One 
can see that for the case of D\ = 0.2D 2 the average DoS 
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in the 7r-band reaches its normal value at fields consider- 
ably smaller than the upper critical field. 

Recently, the c-axis vortex structure in MgB 2 single 
crystals was measured by STM |2^]. Most strikingly, it 
was observed that the spatial extension of the vortex core 
was a few times larger than the length ~10 nm estimated 
from H C 2 . Our model provides an natural explanation of 
this fact. One can see from Fig. |] that the apparent 
vortex size can indeed exceed the size estimated from 
the H C 2 . The magnitude of the enhancement depends on 
the ratio of the diffusion constants in the two bands. As 
follows from numerical calculations, the apparent vortex 
size is approximately given by the expression £ v = 
2.7^2 +0.3£i. The low energy peak (around 2.2 meV) in 
the region between the vortices and its rapid suppression 
by magnetic field is also explained by our model. 

The measured value of £„/£ C 2 = 3 corresponds in our 
model to the ratio = 2. We can not make a 

quantitative comparison between the measured and cal- 
culated values of £„, since scattering parameters in dif- 
ferent bands for MgB2 single crystals of Ref. are not 
known. Moreover, from the available data on resistivity 
and de Haas - van Alphen effect p{| , it follows that the 
cr-band in MgB2 single crystals is in the clean limit. At 
the same time, the 7r-band, probed by c-axis tunneling, 
is moderately dirty, which is consistent with theoretical 
estimates |2l]] . Due to the increase of the effective coher- 
ence length at low energies |2f| , the dirty limit con- 
dition in the 7r-band is certainly satisfied in the energy 
range E < A maX 2. This is consistent with the absence of 
localized states in the vortex core as claimed in Ref. [^3| . 

Our results should still be qualitatively applicable to 
MgB2, even if the cr-band is in the clean limit. Indeed, if 
we focus on the DoS in the 7r-band, which is in the dirty 
limit, then the Usadel equation for Green's function for 
this band is still valid. The only extra input we need is 
the coordinate dependence of the pair potential A2 in the 
7r-band which is coupled to the pair potential Ai in the 
cr-band. The shape of the coordinate dependence of Ai 
does not depend much on the degree of dirtiness, only the 
scale of this dependence changes. Therefore, our result 
can be considered as phenomenological if we define £1 as 
the typical scale of change of Ai. Moreover, £1 has small 
weight in the dependence £1,(^2, 

In conclusion, the vortex core structure in a dirty two- 
band superconductor with weak interband scattering is 
studied theoretically. The distributions of the order pa- 
rameters and local DoS reveal two different spatial scales 
for the two bands, in qualitative agreement with recent 
STM experiments on MgB2. This further supports the 
two band model in MgB2 and also provides an interesting 
new type of vortex core structure. 
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